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As  fuel  cells  become  more  prominent,  new  manufacturing  and  production  methods  will  need  to  be 
developed  to  deal  efficiently  and  effectively  with  increased  demand.  One  necessary  component  of  this 
industrial  growth  is  the  accurate  measurement  of  the  variability  in  the  manufacturing  process.  In  this 
study,  we  present  a  diagnostic  system  that  combines  infrared  thermography  with  a  reactive-flow¬ 
through  technique  to  detect  catalyst-loading  defects  in  fuel-cell  gas-diffusion  electrodes  accurately 
with  high  spatial  and  temporal  resolutions.  Experimental  results  are  compared  with  model  predictions  of 
thermal  response  with  good  agreement.  Data  analysis,  operating-condition  impacts,  and  detection  limits 
are  explored  using  both  experiments  and  simulation.  Overall,  the  results  demonstrate  the  potential  of 
this  technique  to  measure  defects  on  the  millimeter  length  scale  with  temporal  resolutions  appropriate 
for  use  on  a  web-line.  Thus  we  present  the  first  development  stage  of  a  next-generation  non-destructive 
diagnostic  tool,  which  may  be  amenable  to  eventual  use  on  roll-to-roll  manufacturing  lines. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton-exchange-membrane  fuel  cells  (PEMFCs)  are  considered 
a  potential  key  zero-emission  energy-conversion  technology  of  the 
21st  century  [1].  PEMFC  systems  are  currently  being  deployed  in 
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some  early  markets,  and  as  this  commercialization  expands,  there 
must  be  a  transition  from  low-volume,  largely  manual 
manufacturing  methods,  to  high-volume,  continuous  and  auto¬ 
mated  processes.  In  parallel  with  this  transition  to  high-volume¬ 
manufacturing  methods,  there  must  be  an  effort  to  evaluate  and 
develop  quality-control  techniques  that  are  rapid,  noncontact,  and 
non-destructive,  and  which  can  be  implemented  in  continuous  or 
automated  production  lines.  Scalable  quality-control  methods  are 
especially  needed  for  membrane-electrode-assembly  (MEA) 


402 


P.K.  Das  et  al.  /  Journal  of  Power  Sources  261  (2014)  401-411 


components  [2],  including  membranes,  electrodes,  and  gas- 
diffusion  media  (GDM),  which  are  in  many  cases  inspected  visu¬ 
ally  before  use  in  a  stack.  Of  particular  importance,  given  the  high 
cost  of  platinum  and  other  noble-metal  constituents,  is  the  detec¬ 
tion  of  variations  in  catalyst  layers  (CLs).  CL  thickness  variations, 
such  as  a  spot  with  excessive  or  limited  loading,  can  yield  perfor¬ 
mance  losses  beyond  those  expected  from  the  platinum  (Pt)  vari¬ 
ation  or  the  composition  of  the  CL  alone,  especially  at  the  local  area 
near  the  defect  [3-8].  This  could  be  due  to  uneven  compression  of 
the  GDM  at  the  location  of  the  defect,  causing  an  increase  in  serial 
resistance  or  impedance  (e.g.,  contact  losses)  and  leading  to  water- 
management  (e.g.,  flooding  at  the  interface)  and/or  other  issues, 
such  as  cell-performance  reduction  [3,5,9  .  In  addition,  such  vari¬ 
ations  could  lead  to  initiation  sites  for  membrane  and  other  cell 
degradation  and  failure  [10]. 

MEAs  for  PEMFCs  are  typically  fabricated  using  roll-to-roll 
processes,  wherein  a  sheet  of  material  that  is  continuously 
formed  by  coating,  casting,  or  other  methods,  is  conveyed  through 
required  process  steps  by  a  roller  system,  with  the  final  product 
wound  into  a  continuous  roll.  Therefore,  maintaining  the  quality  of 
the  final  product  is  challenging.  For  in-line  inspection,  the  fuel-cell 
industry  is  dependent  on  the  optical  inspection  of  visible  defects  or 
point  measurements  of  coating  thickness.  While  these  measure¬ 
ments  can  be  extremely  useful,  they  will  ultimately  not  be  the  best 
solution  for  continuous  production  of  PEMFC  MEA  components. 
First,  not  all  defects  are  visible,  particularly  relative  to  the  typically 
black,  non-reflecting  electrode  coatings  and  GDM.  Second,  point 
measurements  only  provide  representative  statistical  data  -  they 
cannot  generally  locate  discrete  defects,  even  if  scanned  across  the 
sheet  of  material  being  fabricated,  as  is  often  done  in  the  industry. 
Given  that  discrete  defects  can  act  as  failure-initiation  points  for 
MEAs  in  operation  [11],  it  is  critical  to  develop  an  inspection 
technique  that  measures  every  location  on  the  sheet  material  such 
that  defects  in  MEA  components  can  be  removed  prior  to  assembly 
into  complete  cells  with  minimal  scrap.  This  100%  inspection 
requirement  is  necessary  to  avoid  premature  cell  failures,  which 
often  can  cascade  and  cause  complete  cell  and  stack  failure; 
therefore  necessitating  expensive  replacements.  In  addition,  the 
time  required  for  the  detection  of  defects  is  a  key  metric  for  any 
technique  that  is  required  to  be  on-line  in  a  manufacturing  process 
[2]. 

For  CL  defects,  X-ray  Fluorescence  (XRF)  spectroscopy  has 
recently  been  applied  for  the  in-line  measurement  of  catalyst 
loading  [12,13].  Unfortunately,  this  technique  cannot  at  this  time 
provide  100%  inspection  for  CL  defects.  Acquisition  times  provided 
by  commercially  available  systems  that  provide  quantitative  anal¬ 
ysis  of  Pt  loading  are  such  that,  for  a  roll  width  of  18  in.  and  a  web- 
speed  of  30  ft  min-1,  data  is  only  gathered  on  less  than  1%  of  the 
sample  area.  Therefore,  this  technique  can  only  provide  a  statistical 
metric  of  Pt  loading  for  a  given  production  run.  One  way  around  this 
could  be  a  multi-pass  system  through  the  XRF,  but  this  would  slow 
the  web  speed  and  add  complexity.  Another  option  is  to  use  mul¬ 
tiple  XRF  setups,  but  this  is  cost  prohibitive.  In  summary,  XRF 
spectroscopy  is  an  established  diagnostic  technique,  but  it  has 
limitations  relative  to  use  as  a  diagnostic  tool  for  continuous  pro¬ 
duction  of  PEMFC  electrodes. 

Infrared  (IR)  thermography  may  allow  for  complete,  rapid  (with 
response  times  on  the  order  of  1  s  for  large  areas),  non-contact,  and 
non-destructive  detection  of  defects.  Several  patents  mention  using 
IR  thermography  for  detection  of  various  PEMFC  parameters, 
including  MEA  and  membrane  defects  14-16].  However,  little 
public  information  is  available  on  the  ability  of  specific  excitation 
techniques  that  allow  for  the  detection  of  specific  (size  and  type) 
defects  in  PEMFC  materials.  Recently,  we  have  demonstrated  an  IR 
thermography  and  direct-current  (DC)  electronic-excitation 


method  to  detect  variations  in  CL  thicknesses  on  decals  or  mem¬ 
branes,  i.e.,  electrically  insulating  substrates,  with  high  spatial  and 
temporal  resolution  17].  Although  IR  thermography  techniques 
may  not  be  element  specific  (as  is  XRF),  our  work  shows  that  the  IR- 
DC  technique  can  be  used  to  identify  CL-thickness  variations 
rapidly,  at  a  relatively  low  cost,  while  providing  100%  inspection. 

In  this  paper,  we  examine  an  IR  thermography-reactive  flow 
through  (RFT)  technique  that  is  applicable  to  gas-diffusion  elec¬ 
trodes  (GDEs),  i.e.,  a  CL  coated  on  the  microporous  layer  (MPL)  of  an 
electrically  conductive  gas  diffusion  layer  (GDL).  IR-RFT  detects 
variations  in  Pt  loadings,  which,  in  this  study,  is  identical  to  varia¬ 
tions  in  CL  thickness.  As  an  initial  study  in  the  use  of  IR  detection 
with  reactive  flow-through  for  CL  defects,  this  work  focuses  on 
detecting  square  spots  of  different  thickness  in  a  surrounding  layer 
of  nominal  thickness.  Mathematical  modeling  is  used  to  explore  the 
effects  of  various  process  and  material  parameters. 

The  outline  of  this  paper  is  as  follows.  First,  the  experiment, 
theory,  and  model  for  the  technique  are  described.  Second,  the 
results  of  IR  measurements  with  defined  defects  are  presented  and 
used  for  model  validation.  Third,  both  experiments  and  model  are 
used  to  examine  the  impact  of  different  variables  of  the  technique 
including  gas-mixture  concentration,  flowrate,  defect  size,  and 
aspect  ratio.  Finally,  the  extension  and  applicability  of  this  tech¬ 
nique  on  a  web-line,  including  simulating  a  slot  flow  rather  than  a 
uniform  excitation,  is  discussed. 

2.  Experimental 

IR-RFT  was  employed  to  detect  variations  of  catalyst  (i.e.,  Pt) 
loading  in  GDEs.  The  RFT  technique  involves  the  flow  of  a  dilute 
non-flammable  H2/O2  gas-mixture  through  the  GDE  where  reac¬ 
tion  occurs  at  the  Pt  catalytic  sites  in  the  CL.  For  the  experiments,  a 
flow  of  0.4/0.2%  H2/O2  gas  in  N2  balance  was  introduced  into  an 
inlet  manifold  of  a  bench-top  hardware,  then  passed  through  the 
gas-diffusion  electrode,  and  finally  exited  through  an  outlet 
manifold  as  shown  in  Fig.  1(a).  The  heat  signature  of  the  reaction  at 
the  electrode  was  captured  through  an  IR  transmitting  material 
(POLY  IR®,  Fresnel  Technologies)  using  an  IR  camera  (Jenoptik 
Vario-Cam  HiRes  camera)  with  a  640  x  480  pixel  detector  with  an 
accuracy  of  under  ±1.5  K.  A  uniform  emissivity  of  0.95  was  used  for 
the  GDE.  This  value  was  determined  empirically  by  comparing  a 
piece  of  3  M™  brand  masking  tape,  which  has  a  known  emissivity 
of  0.95  [18],  at  the  same  ambient  conditions  as  the  CL. 

CL  samples  with  and  without  defects  were  prepared  in-house 
from  a  catalyst  ink  containing  20  wt%  Nation  solution  with 
1100  EW  (Ion  Power,  DE2020),  n-propanol,  water,  and  46  wt%  Pt/ 
carbon  catalyst  (TKK,  TEC10E50E).  The  ionomer  to  carbon  mass 
ratio  was  adjusted  to  0.8:1  and  mixing  was  performed  overnight 
using  a  vial,  zirconium  oxide  beads,  and  a  ball  mill.  The  resulting 
catalyst  ink  was  then  sprayed  onto  sheets  of  the  commercially 
available  GDL  material  BC25  from  SGL,  using  a  home-built  ultra¬ 
sonic  spray  system  with  a  Sono-Tek  Impact  spray  head.  The  spray 
system  included  a  syringe  pump,  a  120  kHz  micro-bore  spray  head 
and  power  system,  and  a  programmable  x—y  stage.  The  flowrate  of 
the  catalyst  ink  during  spraying  and  the  speed  of  the  sample  holder 
were  adjusted  to  achieve  a  platinum  loading  of  0.025  mg  cm-2  in  a 
single  coat.  During  spraying,  the  temperature  of  the  sample  holder 
was  controlled  to  80  °C.  No  additional  time  for  drying  was  needed 
between  coats. 

Defects  with  50%  and  100%  loading  reduction  were  intentionally 
introduced  to  select  samples  by  placing  stainless  steel  masks  of 
either  0.25  or  2  cm2  on  the  GDL  material  for  part  of  the  spray 
process.  The  resulting  loading  reduction  of  the  defect  was 
controllable  by  the  number  of  spray  coats  the  mask  was  kept  in 
place.  For  a  loading  reduction  of  50%,  for  example,  the  mask  was 
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Fig.  1.  (a)  Experimental  setup  of  IR-RFT  method  and  (b)  GDE  sample  with  a  square  defect  (bare  spot)  in  the  center  marked  by  dashed  rectangle. 


removed  after  half  the  spray  coats  were  completed,  for  a  loading 
reduction  of  100%  the  mask  was  kept  in  place  during  all  the  spray 
coats.  After  spraying,  the  GDL  sheet  was  cut  into  individual  sample 
pieces,  ensuring  that  all  defects  were  centered  on  the  sample. 

The  bench-top  experimental  hardware  of  the  RFT  technique  and 
an  optical  image  of  a  GDE  with  a  square  defect  in  the  center  are 
shown  in  Fig.  1(a)  and  (b),  respectively.  The  defect  is  highlighted  by 
a  dashed  rectangle  because  the  contrast  between  the  black  CL  and 
the  black  MPL  made  the  defect  difficult  to  see  by  eye. 

3.  Model 

To  guide  the  experiments,  improve  understanding  of  the 
experimental  data,  and  explore  the  technique,  mathematical 
models  of  the  test  setup  were  developed.  Fig.  2(a)  depicts  the 
experimental  setup  and  Fig.  2(b)  shows  the  full  computational 
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domain.  The  computational  domain  is  identical  to  the  electrode 
sample  used  in  the  experiments,  with  the  exception  of  the  brass 
manifolds  shown  in  Fig.  1(a). 


3.1.  Governing  equations 

The  general  energy-balance  equation  used  in  the  model  to  find 
the  temperature,  T,  in  the  porous  GDL,  MPL,  and  CL  is 

(pCp)eq|f  +  pCpU-VT  =  V-(keqVr)  +Q  (1) 

where  p  is  the  gas-mixture  density,  Cp  is  the  gas-mixture  specific 
heat  capacity  at  constant  pressure,  (pCp)eq  is  the  equivalent  volu¬ 
metric  heat  capacity  at  constant  pressure,  T  is  the  temperature,  u  is 
the  superficial  mass-averaged  gas-mixture  velocity,  /<eq  is  the 
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Fig.  2.  (a)  Schematic  of  the  RFT  experimental  setup  with  computational  domain  marked  by  dashed  line,  (b)  Full  computational  domain. 
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equivalent  thermal  conductivity,  and  Q.  is  the  heat  generation  due 
to  the  chemical  reaction  of  O2  and  H2  on  Pt. 

The  equivalent  thermal  conductivity  of  the  porous  matrix/gas 
mixture  system,  /<eq,  is  related  to  the  effective  conductivity  of  the 
porous  layer  (GDL,  MPL,  or  CL),  /<pff,  and  to  the  conductivity  of  the 
gas-mixture,  /q,  by 

ke q  =  (1  -  S)kf  +  skf  (2) 

where  r  is  the  porosity  of  the  porous  layer.  The  effective  thermal 
conductivity  of  the  porous  GDL,  MPL,  and  CL  are  estimated  using  a 
Bruggeman  correction  as  it  provides  a  good  approximation  of  the 
effective  conductivity  [19].  The  equivalent  volumetric  heat  capacity 
of  the  porous  matrix/gas  mixture  system  is  calculated  by 


and  is  nonzero  only  in  the  CL.  The  species  transport  through  the 
porous  GDL,  MPL,  and  CL  is  modeled  using  the  mixture-averaged 
diffusion  model  and  the  species  mass  transport  equation  is  defined  as 

^(pWi)  +  V-(pWjU)  =  -V-;mcU  +  Ki  (10) 

where  Jmd,i  is  the  mass  flux  relative  to  the  mass-averaged  velocity, 
and  R[  is  the  rate  expression  describing  the  production  or  con¬ 
sumption  of  species  ‘i’. 

The  relative  mass  flux  due  to  molecular  diffusion  is  governed  by 
a  Fields  law  expression, 

JnuLi  =  -  (pD^Vw  i  +  p  go™ VM)  (11) 


(pCp)eq  =  (1— £)(pCp)p  +  £(pCp)g  (3) 

where  the  subscripts  p  and  g  represent  the  solid-matrix  of  porous 
layer  and  the  gas-mixture,  respectively.  The  gas-mixture  contains 
hydrogen,  oxygen,  and  nitrogen.  Therefore,  the  gas-mixture  den¬ 
sity,  specific  heat  capacity,  and  thermal  conductivity  are  estimated 
assuming  ideal-gas  properties  by 

p  =  ^Ex M  (4) 

i 

Cp  =  5>iCpi  (5) 

i 

and 

kf  =  (6) 

i 


where  Df1  is  the  mixture-averaged  diffusion  coefficient  and  M  is  the 
mean  molecular  weight  of  the  mixture.  The  molecular  weight  of  the 
mixture  and  the  mixture-averaged  diffusion  coefficient  are  given  by 

,,2) 


and 


Df1 


1  -  oq 


(13) 


respectively,  where  M\  is  molecular  weight  of  species  i,  xi<  the  mole 
fraction  of  species  k,  and  D&  is  the  binary  diffusivity  between 
species  ‘i’  and  ‘k’.  The  rate  expression  for  the  production  or  con¬ 
sumption  of  species  ‘i’  is  defined  as 

Ki  =  Mjfir  (14) 


respectively,  where  R  is  the  ideal-gas  constant,  and  Xi,  oq,  Mi,  and  k\ 
are  the  mole  fraction,  mass  fraction,  molecular  weight,  and  thermal 
conductivity  of  species  ‘i\  respectively.  It  is  worth  to  note  that  the 
gas-mixture’s  specific  heat  capacity  is  estimated  using  the  mass- 
weighted-mixing-law  and  the  thermal  conductivity  is  estimated 
using  a  simplified  molar-weighting  approach. 

The  heat  generation  term  is  described  by  the  reaction  of 
hydrogen  as 

Q.  =  Krfrcomb  (?) 

where  Rr  is  reaction  rate  and  /iCOmb  is  the  lower  heating  value  of 
hydrogen  (119.96  MJ  kg-1).  The  reaction  rate  is  described  with  the 
one-step  mechanisms  of  Mitani  and  Williams  [20]: 

Rr  =  AvBCH2  exp(-Tj  (8) 

with  B  =  1.4  x  103  cm  s-1  and  E  =  14.9  x  103  J  mol-1.  Here  Av  is  a 
fitting  parameter  to  account  for  interfacial  surface  area  to  volume 
ratio. 

The  gas-mixture  transport  through  the  porous  GDL,  MPL,  and  CL 
is  subject  to  the  following  mass-conservation  equation,  which  in¬ 
corporates  Darcy’s  law  for  flow  in  the  porous  media, 


3.2.  Boundary  conditions  and  solution  methodology 

The  above  governing  equations  require  boundary  conditions. 
Since  the  brass  manifolds  are  not  explicitly  modeled,  the  heat-sink 
effect  of  the  brass  manifolds  was  accounted  for  using  an  isothermal 
boundary  condition,  which  was  found  to  have  a  negligible  impact 
on  the  numerical  temperature  profile  because  the  defects  studied 
are  located  at  the  center  of  the  sample.  Between  solid  and  gas 
phases,  the  model  assumes  local  thermal  equilibrium. 

For  the  energy  equation,  a  thermal  boundary  condition  at  the  inlet 
manifold/GDL  interface  is  specified  by  the  gas-mixture  temperature, 

T  =  Tin  for  6£ez  =  0  (15) 

where  9Q  represents  the  boundary  of  the  computational  domain.  At 
the  outlet  manifold/CL  interface,  a  convective  cooling  boundary 
condition  is  assigned, 

n-(-kaVT)  =  hConv(T  —  TYoom)  for  dQ&z  =  ti  (16) 

where  ti  is  the  combined  thickness  of  GDL,  MPL,  and  CL  layer  and 
hc onv  is  the  equivalent  convection  heat-transfer  coefficient  that  is 
estimated  using  an  expression  similar  to  Handley  and  Heggs 
expression  21]: 


a^ 

at 


(£/>)  +  V- 


—  Qm 


(9) 


% 


0.255/<f 


Pr^Re^ 


(17) 


where  p,  /i,  k,  and  p  are  the  density,  dynamic  viscosity,  permeability, 
and  pressure,  respectively.  The  right  side  of  above  equation  is  the 
source  term  that  accounts  for  mass  change  due  to  the  cell  reaction 


where  Pr  and  Re d  are  the  Prandtl  number  of  the  gas-mixture  and 
Reynolds  number  based  on  manifold  inlet  diameter,  respectively, 
and  af  is  a  fitting  parameter  which  is  fit  once  and  then  not  varied 
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Fig.  3.  Steady  state  heat  signatures  at  flowrates  of  (a)  500  seem,  (b)  1000  seem,  and  (c)  1500  seem.  The  GDE  used  is  pristine  with  a  loading  of  0.1  mg-Pt  cm-2  and  0.4%  H2  and  0.2%  02 
in  N2  gas-mixture.  The  dashed  rectangle  in  part  (a)  is  the  data  collection  window,  which  is  about  75%  of  the  total  sample  area. 


(a)  500  seem  (b)  1 000  seem  (c)  1 500  seem 


among  the  different  simulation  cases.  The  Prandtl  number  of  the 
gas-mixture  is  estimated  using  the  mixture  composition,  while  the 
Reynolds  number  is  estimated  using  the  following  expression. 


Reo 


QmD 

~W 


(18) 


where  qm  is  the  gas  flowrate,  D  is  the  inlet-tube  diameter,  A  is  the 
area,  and  v  is  the  kinematic  viscosity. 

For  the  mass  conservation  equation,  inlet  mass-flux  is  defined 
based  on  the  gas  flowrate,  while  a  pressure  boundary  condition  is 
applied  the  outlet  manifold/CL  interface  as 

p  =  Proom  for  a£ez  =  tL  (19) 

For  the  species  transport  equation,  inlet  mass-fractions  are 
defined  based  on  the  gas-mixture  composition,  while  a  Neumann 
boundary  condition  is  applied  at  the  outlet  manifold/CL  interface 

V^i  =  0  for  0 Q&z  =  tL  (20) 

For  all  other  boundaries  an  insulation  or  a  symmetry  boundary 
condition  is  applied. 

The  governing  equations  were  solved  using  commercial  finite- 
element  software  COMSOL  Multiphysics®.  The  3D  computational 
domain  was  discretized  into  a  tetrahedral  mesh  and  second-order 
Lagrangian  elements  were  used  with  the  Direct  solver.  The  solu¬ 
tions  were  considered  converged  when  the  tolerance  value  was 
below  1CT5  and  were  mesh  independent  (at  least  30  x  35  x  35 
elements). 


4.  Results  and  discussion 

4.1.  Pristine  samples  and  model  validation 

To  understand  the  viability  of  the  technique,  pristine  GDEs,  i.e., 
GDEs  with  nominally  uniform  Pt  loading  were  fabricated  and  tested 
in  the  setup  shown  in  Fig.  1.  Fig.  3  shows  IR  thermographs  of  a 
pristine  sample  at  steady  state  over  a  range  of  flowrates.  It  is 
apparent  that  the  heat  signature  from  the  CL  is  readily  detectable, 
owing  to  the  large  heat  of  reaction  of  hydrogen.  It  should  be  noted 
that  IR  thermography  is  essentially  a  surface  technique,  although 
the  thinness  of  the  CL  ensures  that  the  temperature  is  uniform 
through  the  thickness  of  the  layer.  The  heat  of  reaction  is  removed 
through  the  top  and  bottom  CL  surfaces,  by  convection  with  the  gas 
stream,  and  by  conduction  to  the  brass  manifolds.  This  resulting 
thermal  balance  results  in  Fig.  3.  The  in-plane  temperature  gradient 
observed  in  the  thermographs  is  due  mainly  to  the  presence  of  the 


large  thermal  mass  of  the  brass  manifold  that  essentially  induces  an 
isothermal  boundary,  and  may  also  be  due  to  some  channeling  of 
the  flow  in  the  middle  of  the  domain.  Fig.  3  also  shows  the  increase 
in  heat  signature  from  31  to  36  °C  as  the  gas  mixture  flowrate  is 
increased  three-fold.  The  increased  flowrate  corresponds  to  an 
increased  flux  and  thus  rate  of  reaction  of  the  hydrogen. 

To  understand  the  impact  of  flowrate  in  more  detail,  Fig.  4 
shows  the  maximum  and  average  temperature  rise  for  the  pris¬ 
tine  sample  as  a  function  of  flowrate.  The  maximum  temperature 
represents  the  highest  temperature  measured,  while  the  average 
temperature  represents  the  mean  temperature  within  a  collection 
window  of  about  75%  of  the  total  sample  area  around  the  center  of 
the  sample  as  shown  in  Fig.  3(a).  The  size  of  the  collection  window 
was  kept  smaller  than  the  sample  area  to  minimize  the  edge  effects 
observed  in  Fig.  3.  The  two  temperatures  track  consistently,  which 
is  expected  since  the  data  originate  from  the  same  set  of  experi¬ 
ments.  The  shape  of  the  curves  in  Fig.  4  is  the  result  of  two 
competing  processes.  Below  1500  seem  the  temperature  increases 
with  increasing  flowrates  due  to  an  increase  in  hydrogen  flux  and 
thus  a  higher  reaction  rate.  However,  as  the  flowrate  becomes 
larger  than  1500  seem  a  plateau  is  reached.  This  was  attributed  to 
complete  utilization  of  all  active  Pt  catalyst  sites  in  the  CL.  In  other 
words,  at  gas  flows  above  1500  seem  hydrogen  is  flowing  through 
the  CL  without  reacting.  At  flowrates  above  2500  seem  the 


Fig.  4.  Measured  temperature  rise  as  a  function  of  gas  flowrate  for  pristine  GDE  with 
Pt  loading  of  0.1  mg-Pt  cm-2  and  0.4%  H2  and  0.2%  02  in  N2  gas-mixture.  Error  bars 
show  90%  confidence  intervals  for  seven  total  runs. 
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temperature  decreases  due  to  what  we  believe  to  be  increased 
convective  cooling  of  the  GDE. 

To  explore  the  energy  balance  and  parameter  space  in  more 
detail,  a  mathematical  model  was  formulated  as  described  above. 
The  model  was  validated  by  comparison  of  the  temperature  rise  of 
different  pristine  GDE  samples  as  shown  in  Fig.  5  for  various 
nominal  Pt  loadings  and  hydrogen  concentrations  at  a  flowrate  of 
1500  seem.  As  shown,  the  model  predictions  demonstrate  good 
agreement  with  the  experimental  data  for  various  Pt  loadings  with 
a  few  percent  error.  The  data  in  Fig.  5(a)  indicate  increased  extent  of 
reaction  of  hydrogen  with  increasing  Pt  loading  and  confirm  that  at 
these  specific  experimental  conditions,  and  over  the  range  of 
loadings  studied,  not  all  hydrogen  is  reacting  and  contributing  to 
the  heat  signature  of  the  GDE. 

Fig.  5(b)  shows  the  impact  of  hydrogen  concentration  in  the 
reactive  gas  mixture  on  the  temperature  response.  For  this  effect, 
the  value  comparison  between  model  and  experiment  are  in  even 
closer  agreement  than  the  comparison  in  Fig.  5(a),  which  also 
highlights  the  fact  that  the  model  does  not  require  re-calibration 
for  the  different  concentrations.  It  is  also  apparent  that  the  tem¬ 
perature  response  is  highly  sensitive  to  the  hydrogen  concentra¬ 
tion,  where  going  towards  2%  provides  over  a  six-fold  increase  in 
the  temperature  rise;  it  should  be  noted  that  2%  is  still  only  half  of 
the  flammability  limit  of  hydrogen  in  air. 

4.2.  Detection  of  defects 

In  this  section  we  examine  the  feasibility  of  the  introduced  IR- 
RFT  method  to  detect  defects  that  were  intentionally  introduced 
into  the  CL.  Square  defects  of  different  sizes  and  degrees  were 
tested  and  square  and  rectangular  defects  simulated  with  different 
flowrates,  nominal  Pt  loadings,  and  hydrogen  concentrations. 

Fig.  6(a)— (c)  shows  IR  thermographs  and  corresponding  tem¬ 
perature  profiles  of  5  cm  x  5  cm  GDE  samples  with  square  defects 
at  the  center  of  the  sample.  Data  is  shown  for  a  gas  flowrate  of 
1500  seem  at  10  and  20  s  as  thermographs  in  the  left  and  center, 
respectively,  and  as  temperature  profiles  along  the  dashed  line  in 
the  right.  The  nominal  Pt  loading  of  the  pristine  area  was 
0.2  mg  cm-2.  Fig.  6(a)  shows  results  for  a  2  cm2  square  100%  loading 
reduction  defect,  i.e.  a  2  cm2  area  without  any  Pt  loading.  The  defect 
was  indicated  after  10  s  by  a  temperature  reduction  of  over  4  °C, 
which  is  apparent  in  the  thermographic  image  on  the  left  as  well  as 
the  temperature  profile  on  the  right. 

In  Fig.  6(b),  results  are  shown  for  a  0.25  cm2  defect  with  a  100% 
loading  reduction.  The  temperature  variation  at  the  defect  location 


is  about  1  °C.  Thus,  the  defect  is  detectable  but  it  could  be  possibly 
enhanced  by,  for  example,  increasing  the  gas  flowrate  as  indicated 
by  Fig.  4  or  by  increasing  the  hydrogen  concentration  as  indicated 
by  Fig.  5(b).  A  smaller  defect,  in  this  case  a  50%  loading  reduction  of 
the  same  size,  becomes  more  challenging  to  detect  as  shown  in 
Fig.  6(c).  Although  after  20  s  the  defect  signature  is  visible  in  the 
temperature  profile,  the  temperature  differential  is  very  small. 

To  enhance  the  signal,  as  mentioned,  a  higher  flowrate  of  inlet 
gas  could  potentially  be  used.  Fig.  7(a)  and  (b)  shows  the  maximum, 
minimum,  and  differential  temperature-time  profiles  of  a 
5  cm  x  5  cm  defect  GDE  sample  with  a  1  cm2  50%  defect  at  the 
center  for  gas  flows  of  1500  and  5000  seem,  respectively.  The 
temperature  rise  for  5000  seem  is  faster  than  for  the  1500  seem 
case,  thus  showing  that  one  can  obtain  response  times  on  the  order 
of  a  few  seconds.  However,  the  temperature  increase  is  larger  for 
the  1500  seem  case  because  of  the  increased  convective  cooling  at 
5000  seem. 

As  seen  in  Fig.  7,  one  must  be  cognizant  of  choice  of  design 
criteria  since  just  using  the  maximum  temperature  is  not  neces¬ 
sarily  the  best  metric  since  variations  in  defects  and  loadings  could 
give  profiles  that  fluctuate  and  thus  a  threshold  criterion  may 
become  ill-defined.  Similarly,  the  temperature  difference,  if  used, 
must  be  carefully  defined  with  respect  to  a  reference  or  average 
temperature. 

To  understand  the  impact  of  flowrate  in  more  detail,  simulations 
at  different  flowrates  were  conducted.  Fig.  8  shows  the  steady-state 
temperature  profiles,  the  extent  of  reaction  for  the  ^/C^mixture, 
and  the  pressure  drop  through  the  GDE.  The  steady-state  temper¬ 
ature  profiles  are  shown  from  the  middle  of  one  edge  to  the  center 
of  a  5  cm  x  5  cm  electrode  with  a  1  cm2  square  defect  at  the  center 
for  various  gas-mixture  flowrates  and  for  Pt  loading  of  (a)  0.2  mg- 
Pt  cm-2  and  (b)  0.4  mg-Pt  cm-2.  The  extents  of  reaction  are  shown 
as  the  percentage  hydrogen  mass  fractions  at  the  exit  of  the  GDE  for 
various  flowrates  and  for  Pt  loading  of  (c)  0.2  mg-Pt  cm-2  and  (d) 
0.4  mg-Pt  cm-2.  The  pressure  drop  (difference  between  the  local 
pressure  and  ambient  pressure)  data  in  Fig.  8(e)  are  presented 
along  the  thickness  of  the  sample  for  two  flowrates  under  the 
defect  area  and  the  pristine  area  for  Pt  loading  of  0.2  mg-Pt  cm-2. 
Similar  to  the  experimental  data,  the  modeling  results  show  that 
the  temperature  in  the  pristine  area  is  always  higher  than  the 
temperature  in  the  defect  area.  Also,  a  higher  loading  always  gen¬ 
erates  a  higher  temperature  due  to  more  complete  combustion  of 
the  hydrogen.  As  shown  in  Fig.  8(c)  and  (d),  the  exit  mass-fraction 
percentage  of  hydrogen  is  significantly  higher  for  lower  loadings.  In 
fact,  doubling  of  the  loading  doubled  the  differential  temperature. 


Fig.  5.  Comparisons  between  experimental  measurements  and  model  predictions  for  5  cm  x  5  cm  pristine  samples  at  a  gas  flowrate  of  1500  seem:  Temperature  rise  as  a  function  of 
(a)  nominal  loading  for  0.4%  H2  and  0.2%  02  in  N2,  and  (b)  hydrogen  concentration  for  a  nominal  Pt  loading  of  0.4  mg-Pt  cm-2' 
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Fig.  6.  IR  thermographs  at  10  and  20  s,  and  corresponding  line  temperature  profiles  of  5  cm  x  5  cm  defect  samples  with  0.2  mg-Pt  cm  2  nominal  loading  at  1500  seem  of  0.4%  H2 
and  0.2%  02  in  N2:  (a)  2  cm2  defect  of  100%  loading  reduction  (i.e.  a  bare  spot),  (b)  0.25  cm2  defect  of  100%  loading  reduction,  and  (c)  0.25  cm2  defect  of  50%  loading  reduction. 


Fig.  7.  Temperature  profiles  of  a  5  cm  x  5  cm  GDE  sample  with  0.2  mg-Pt  cm  2  nominal  loading  and  1  cm2  50%  defect  using  0.4%  H2  and  0.2%  02  in  N2  at  (a)  1500  and  (b)  5000  seem, 
respectively. 
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Fig.  8.  Simulated  (a,  b)  temperature  profiles,  (c,  d)  extents  of  reaction,  and  (e)  pressure  drop  through  a  5  cm  x  5  cm  GDE  with  1  cm2  50%  defect  with  various  gas-mixture  flowrates 
using  0.4%  H2  and  0.2%  02  in  N2.  The  temperature  profiles  are  from  CL  edge  to  defect  center  for  loadings  of  (a)  0.2  mg-Pt  cm-2  and  (b)  0.4  mg-Pt  cm-2,  the  extents  of  reaction  are 
presented  as  the  exit  mass-fraction  percentage  of  hydrogen  for  (c)  0.2  mg-Pt  cm-2  and  (d)  0.4  mg-Pt  cm-2,  and  the  pressure-drop  data  are  for  0.2  mg-Pt  cm-2. 


Fig.  8  reveals  that  higher  gas  flowrates  result  in  lower  temperatures 
in  both  the  pristine  and  defect  areas.  However,  the  differential 
temperatures  between  the  pristine  and  defect  areas  seem  almost 
independent  of  the  gas  flowrate,  which  is  a  key  finding  in  terms  of 
optimizing  detection  and  response  time.  In  other  words,  high 
flowrates  enable  faster  response  times  because  reactants  reach  the 
catalyst  layer  faster  (see  Fig.  7),  while  still  maintaining  a  detectable 
differential  temperature.  The  pressure-drop  data  reveal  that  the 


entire  pressure  drop  is  happening  between  the  MPL  and  CL  with  a 
negligible  pressure  drop  in  the  GDL.  The  simulation  data  also  reveal 
that  heat  convection  and  conduction  are  not  dominant  in  the 
experimental  set-up  as  the  temperature  is  within  1  °C  of  the  pre¬ 
dicted  adiabatic  temperature  rise  (i.e.,  assuming  no  heat  transfer 
out  of  the  catalyst  layer).  For  instance,  the  adiabatic  temperature 
rise  in  the  pristine  are  for  0.2  mg-Pt  cm-2  loading  and  1000  seem  is 
about  6.5  °C,  whereas  the  actual  temperature  rise  is  5.7  °C. 
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However,  the  variation  between  the  adiabatic  and  the  actual  tem¬ 
perature  rise  is  slightly  higher  at  higher  flowrates  due  to  higher 
convective  heat  removal. 

Using  the  model,  one  can  predict  the  detection  limit  of  the 
technique  relative  to  different  degrees  of  defect,  defined  as  the 
percent  of  loading  reduction,  where  0%  loading  reduction  indicates 
a  pristine  electrode  and  100%  loading  reduction  indicates  a  bare 
spot.  Fig.  9  shows  the  differential  temperature  between  the  pristine 
and  defect  areas  as  a  function  of  defect  degree.  Interestingly,  the 
data  is  almost  linear  with  defect  degree.  The  lower  flowrate  pro¬ 
vides  a  slightly  higher  temperature  rise  than  the  higher  flowrate, 
which  is  due  to  the  convective  cooling  effect  at  the  higher  flowrate 
(see  Fig.  4).  However,  when  using  the  detection  limit  we  defined  in 
our  previous  work  for  IR-DC  [17],  i.e.  AT=  1  °C,  the  effect  of  flowrate 
is  minor.  The  effect  of  flowrate,  and  thus  convective  cooling,  be¬ 
comes  more  prominent  as  reduction  of  loading  in  the  defect  area 
increases.  Fig.  9  also  reveals  that,  for  this  defect  and  with  a  1  °C 
detection  limit,  the  limitation  on  detecting  the  degree  of  the  defect 
is  when  the  loading  reduction  in  the  defect  area  is  less  than  22%  of 
the  nominal  loading. 

The  data  shown  in  Fig.  5  predicted  that  the  thermal  response  of 
a  pristine  GDE  can  potentially  be  increased  by  increasing  the 
hydrogen  concentration.  Fig.  10  shows  a  similar  prediction  for  the 
differential  temperature  caused  by  a  defect.  The  simulations  indi¬ 
cate  that  by  increasing  the  hydrogen  concentration  in  the  inlet  gas 
stream,  smaller  defects  may  be  able  to  be  detected.  Clearly,  higher 
hydrogen  concentrations  raise  the  reaction  rate  and  generate  a 
greater  differential  temperature.  Therefore,  a  defect  that  is  not 
detectable  with  0.4%  hydrogen  (see  Fig.  9)  may  be  able  to  be 
detected  by  increasing  hydrogen  concentration  in  the  inlet  gas- 
stream  (of  course  for  safety  reasons  one  wants  to  remain  below 
the  flammability  limit  of  4%  H2).  In  fact,  increasing  the  hydrogen 
concentration  is  preferable  to  increasing  the  hydrogen  gas  flow, 
because  the  latter  increases  the  amount  of  counterproductive 
convective  cooling.  Increasing  the  hydrogen  concentration  five¬ 
fold,  from  0.4  to  2%,  provides  an  almost  seven-fold  increase  in 
temperature  difference  for  this  defect  at  a  flowrate  of  1000  seem. 
Thus,  when  using  2%  H2  and  1%  O2  in  N2,  the  sensitivity  for 
detecting  defects  will  be  enhanced  and,  assuming  a  detection 
criteria  of  1  °C  variation,  smaller  defects  may  become  detectable 
(e.g.,  simulations  show  down  to  0.0625  cm2  for  bare  spots). 

To  understand  the  impact  of  defect  shape  (such  as  non-square 
defect  or  micro-crack),  simulation  results  are  presented  in  Fig.  11 
that  study  the  effect  of  the  defect  aspect  ratio  (AR),  which  is 


Fig.  9.  Temperature  rise  between  the  defect  and  non-defect  area  as  a  function  of  defect 
loading  at  two  different  flowrates  for  a  5  cm  x  5  cm  GDE  sample  with  a  1  cm2  defect. 


Fig.  10.  Temperature  differential  between  the  defect  and  non-defect  (pristine)  area  as 
a  function  of  hydrogen  concentration  for  a  5  cm  x  5  cm  GDE  sample  with  a  1  cm2  50% 
loading  reduction  defect  at  two  different  flowrates. 

defined  as  the  height  (H)  to  width  (W)  ratio  of  the  defect.  Fig.  11(a) 
shows  the  temperature  profiles  across  the  electrode  along  the 
white  dashed  line  shown  in  the  simulated  thermographic  image 
shown  in  Fig.  11(b)  and  (c)  shows  the  temperature  differential  be¬ 
tween  defect  and  non-defect  area  for  a  fixed  defect  height  of  5  mm 
and  different  defect  widths  ranging  from  5  mm  to  200  pm.  As 
observed  in  the  temperature  profiles,  the  maximum  temperature  in 
the  pristine  area  remains  identical  for  all  the  aspect  ratios;  how¬ 
ever,  the  temperature  in  the  defect  region  decreases  as  the  aspect 
ratio  increases.  The  results  of  Fig.  11(c)  show  that  for  narrower 
defect  geometries  such  as,  for  example,  a  micro-crack 
(200  pm  x  5  mm  crack,  AR  =  25),  the  temperature  rise  between 
the  defect  and  non-defect  area  decreases.  It  is  further  predicted  that 
for  such  a  micro-crack,  the  temperature  differential  is  larger  than 
3  °C  when  using  2%  H2  containing  gas.  Therefore,  the  current 
technique  may  be  capable  of  detecting  micro-cracks  when  using 
higher  hydrogen  concentrations.  In  theory,  the  IR-detection  sensi¬ 
tivity  can  be  as  good  as  a  few  micrometers  and  with  use  of  the 
higher  concentrations  and  optimal  flowrate,  the  micro-cracks  and 
micro-defects  that  are  expected  to  be  troublesome  in  the  fuel  cells 
may  possibly  be  detected. 

5.  Implementation  on  web-line 

5.1.  Flow  arrangement 

The  experimental  hardware  introduced  in  this  work  is  the  first 
experimental  development  step  towards  an  IR-RFT  based  diag¬ 
nostic  system  that  is  envisioned  to  operate  on  an  industrial  pro¬ 
cessing  line.  One  common  method  to  apply  gases  in  such  a  web 
environment  is  through  the  use  of  a  gas  knife,  which  transports  the 
gas  through  either  a  slit  or  an  array  of  holes  onto  the  moving  web  to 
create  an  impinging  flow.  Thus,  instead  of  a  uniform  flux  as  shown 
in  Fig.  1  throughout  the  domain  of  interest,  there  would  be  a  more 
complicated  flow  distribution  onto  and  around  the  GDE.  While  an 
analysis  of  this  situation  is  beyond  the  scope  of  this  paper,  it  is 
worthwhile  to  examine  if  such  a  slit  or  line  flow  through  the  GDE  is 
feasible  for  the  detection  of  defects.  The  setup  in  Fig.  1  and  the 
model  domain  in  Fig.  2  were  modified  as  shown  in  Fig.  12(a)  such 
that  flow  only  occurred  through  a  slit.  The  slit  was  created  by 
adding  a  PTFE  sheet  that  contained  a  3  mm  wide  rectangular  cut¬ 
out  to  the  inlet-side  of  the  GDE.  From  the  thermographic  image 
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Fig.  11.  Modeling  predictions  for  a  gas  mixture  of  2%  H2  and  1%  02  in  N2,  a  variety  of 
geometric  aspect  rations  and  flows,  (a)  Effect  of  geometric  aspect  ratio  of  the  defect  on 
the  temperature  profiles  across  the  electrode  along  the  white  dashed  line  for  a  50% 
loading  and  0.5  cm  wide  defect,  (b)  Modeling  results  of  thermographic  image  with 
geometric  aspect  ratio  of  1.  (c)  Temperature  differential  between  defect  and  non-defect 
area  as  a  function  of  aspect  ratio  at  two  different  flowrates. 


shown  in  Fig.  12(b),  which  shows  the  results  of  a  2  cm2  100% 
reduction  defect,  it  is  apparent  that  the  bare  spot  is  still  detectable. 
Thus,  the  analyses,  results,  and  discussion  presented  above  for  a 
flow-through  technique  leads  us  to  believe  that  a  suitably  modified 


technique  is  promising  for  implementation  on  a  web-line.  The 
impact  of  impinging  flow  and  flow  penetration  into  the  CL,  rather 
than  through  the  GDL  and  MPL,  will  be  critical  to  a  web-line 
implementation  and  is  currently  being  investigated. 

5.2.  Non-destructiveness 

Temperature  rises  have  been  observed  and  predicted  that 
ranged  up  to  35  °C  for  0.4%  H2  and  to  70  °C  for  2%  H2,  respec¬ 
tively.  While  these  temperature  increases  specifically  for  the 
latter  concentration  may  appear  large,  they  are  still  below  typical 
operating  temperatures  of  PEMFCs.  For  example,  the  most  sus¬ 
ceptible  component  may  be  the  ionomer,  but  the  a-transition 
temperature  of  Nation  is  over  100  °C  [22]  and  thus  morphological 
changes  are  not  expected  even  for  the  case  with  2%  or  slightly 
higher  hydrogen  concentrations.  In  addition,  both  the  experi¬ 
mental  and  numerical  results  shown  are  for  non-moving  mate¬ 
rials  that  are  continuously  excited,  i.e.  for  up  to  20  s  in  the 
experimental  data  and  steady  state  solutions  for  the  numerical 
models.  Application  of  this  technique  in  a  web-line  environment 
involves  diagnosing  moving  GDE  substrates.  Thus  excitation  of 
any  point  of  the  GDE  is  only  applied  for  a  few  seconds  and  the 
temperature  rise  is  expected  to  be  lower  than  reported  here.  The 
final  temperature  rise  on  a  web-line  will  be  a  function  of  web 
speed,  gas  flow,  gas  composition,  and  potentially  other  parame¬ 
ters.  Thus,  the  technique  is  expected  to  be  non-destructive  to  the 
catalyst  and  the  GDE  in  general. 


6.  Conclusions 

Quality-control  diagnostics  of  fuel-cell  materials  and  compo¬ 
nents  are  required  for  high-volume  production  of  fuel  cells.  The 
first  development  stage  of  a  non-destructive,  rapid  technique  to 
measure  Pt-loading  variations  in  mass-produced  gas-diffusion 
electrodes  was  introduced.  This  development  stage  of  the  diag¬ 
nostic  technique  involves  flowing  a  dilute  hydrogen/oxygen  gas 
mixture  through  the  gas-diffusion  electrode,  and  then  measuring 
the  subsequent  temperature  response  due  to  the  reaction  of  the 
hydrogen  with  oxygen  on  Pt  using  infrared  thermography.  Both 
experiments  and  modeling  demonstrate  that  this  technique  allows 
one  to  resolve  spatially  the  variations  in  the  catalyst-layer  loadings 
by  thermographic  imaging,  where  the  temperature  response  is  a 
function  of  the  Pt  loading. 

The  technique  was  successfully  employed  to  detect  a  variety  of 
catalyst  layer  defects  that  were  intentionally  introduced  into  the 
GDE.  The  signal  strength  and  speed  for  defect  recognition  was 
found  to  be  dependent  on  the  applied  operating  parameters. 
Increasing  the  inlet  gas  flowrate  reduces  the  time  required  for 
detecting  defects,  however,  increased  flowrates  were  also  observed 
to  result  in  convective  cooling  and  a  decrease  in  the  maximum 
temperature.  The  hydrogen  concentration  of  the  gas  mixture  was 
identified  as  a  key  parameter  for  using  the  technique.  Modeling 
work  indicated  that  for  the  given  operating  conditions  the  tem¬ 
perature  response  could  be  enhanced  by  10  °C  or  more  when 
switching  from  a  hydrogen  concentration  of  0.4%  to  concentrations 
of  2%  or  higher. 

The  results  demonstrated  on  the  first  development  stage  of  this 
technique  show  great  promise  for  further  development  to  an  in¬ 
line  quality-control  diagnostic  for  fuel-cell  gas-diffusion  elec¬ 
trodes.  The  technique  is  cost  effective,  rapid,  allows  for  100%  areal 
inspection,  and  utilizes  only  non-hazardous  gas,  i.e.,  a  H2/02  gas 
composition  that  is  well  below  the  flammability  limit.  Future  work 
will  focus  on  the  second  development  stage  of  this  technique:  the 
transition  from  a  static  sample  chamber  with  non-moving  sample 
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Fig.  12.  (a)  Experimental  configuration  for  line  flow  experiment  using  a  slotted  gasket,  (b)  Thermographic  image  of  a  2  cm2, 100%  reduction,  0.2  mg-Pt  cm  2  nominal  loading  sample 
at  a  flowrate  of  3000  seem  and  hydrogen  concentration  of  0.4%. 


materials  shown  in  this  work  to  moving  sample  materials  in  an 
open  exposed  environment. 
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